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3-D KINEMATICS AND TOLERANCE VARIATION ANALYSIS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is related to co-pending U. S. Provisional Application for Patent Serial No . 
60/220,897, entitled: SIGMUND 3D - KINEMATICS MODULE ALGORITHM, and filed July 
26, 2000, which is hereby incorporated by reference as though set forth in full herein. 

TECHNICAL FIELD OF THE INVENTION 

The principles of the present invention generally relate to a kinematic modeling system, and 
more specifically, but not by way of limitation, to a dynamic modeling system including kinematic 
modeling having tolerances. 

BACKGROUND OF THE PRESENT INVENTION 

Different approaches for kinematics analysis have been used in different mechanical design 
applications. Approaches for modeling an open-loop system are much simpler and more 
straightforwardthanmodeling a closed-loop system. In modeling an open-loop system, generally, 
there is no need to solve non-linear system simultaneous equations in such a system. For a closed- 
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loop system, however, solving non-linear system simultaneous equations is required, and, thus, an 
eflBcient and robust algorithm becomes important. One popular approach for modeling a closed-loop 
system is using a vector loop approach, which requires forming at least one closed-loop in a 
mechanical system. However, the vector loop approach cannot be applied without knowing whether 
5 an open-loop or closed-loop system exists. 

In performing the vector loop approach, kinematic governing equations are obtained by 
forming the proper number of closed-loops in a mechanical system. In general, a vector equation 
(which is equivalent to three scale equations) is obtained for each closed-loop under consideration. 
Determining the number of closed-loops and which loop is required for obtaining solutions depends 

10 on the configuration of the mechanical system, i.e., how joints/drivers are arranged in a loop. 
However, if a mechanism becomes complicated, the kinematic governing equations v^U typically 
become difficult to obtain. Thus, using the vector loop approach is not particularly efficient or robust. 
Furthermore, whenintegratmgthe vector loop approach with variation analysis, the solution will be, 
undoubtably, difficult to converge (especially with large variations) because, generally, there are no 

1 5 means available to control potential axis misalignments in the mechanisms due to variations. One 
typical example of a mechanism is a planar 4-bar linkage. It has one degree of freedom in a plane, 
but has zero degree of fi^eedom in space. Consequently, if variations cause the linkages to he on 
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diflFerent planes, the solution will almost never converge to areasonable accuracy. For this reason, 
the current available vector loop approach, adopted in kinematics analysis with variation analysis 
packages, cannot actually solve the kinematic governing equations with variation conditions. 
Alternatively, it is a common to practice to estimate possible results based on kinematics constraints, 
5 Therefore, adopting another approach is required to solve the kinematic governing equations with 
variation conditions. 

On the other hand, in perfbrmingtheEuler parameters approach, the governing equations are 
obtained by equating joint constraints directly through a coordinate transformation using a set of 
nonlinear equations expressed by general coordinates of two connected parts. Adopting Euler 

1 0 parameters in a transformation matrix simplifies the kinematic governing equations, and also takes into 
account the flexibility of a piece part. Because the kinematic governing equations are obtained by 
equating the joint constraints directly through the coordinate transformation, potential misalignments 
in the mechanisms due to variations can be easily controlled . Moreover, because there is generally 
no need to form any closed-loops, the Euler parameters approach can thus be appUed in both open- 

1 5 loop and closed-loop systems without having to know the type of the system. Additionally, the 
kinematic governing equations obtained are independent of the configuration of the mechanical 
system, i.e., howjoints/drivers are arranged in the system. Therefore, even though a mechanism 
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becomes very complicated (i.e., having multiple loops and/or multiple joints in a piece part and/or 
multiple piece parts connected to a joint), the kinematic governing equations can still be easily 
obtained. Thus, it can be said that the Euler parameters approach is efficient and robust, and also 
proper for being integrated with variation analysis packages. 
5 Presently, conventional modeling systems generate kinematic models of a dynamic mechanical 

system and apply random tolerance parameters at different steps along atrajectory to determine the 
results of the modeled mechanical system. By applying random tolerance parameters at different 
steps, the piece parts are changed throughout the trajectory and a user is unable to interactively view 
and/or obtain results for a single mechanical system. 



10 SUMMARY OF THE INVENTION 

To overcome the problems of conventional modeling techniques, integration of kinematic and 
tolerance modeling is performed to allow interactive viewing and/or obtaining results of a dynamic 
mechanical system. One embodiment of the present invention includes a system and method for 
providing graphical and analytical dynamic modeling of a multi-dimensional mechanical mechanism. 
1 5 The system includes at least one processor for computing dynamic simulation of movement of the 
multi-dimensional mechanical mechanism, and rendering the multi-dimensional mechanical mechanism. 
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The dynamic simulation simultaneously performs kinematic and tolerance modeling. The system 
further includes a display for displaying the rendered multi-dimensional mechanical mechanism during 
the dynamic simulation. 



BRIEF DESCRIPTION OF THE DRAWINGS 

5 A more complete understanding of the method and apparatus of the present invention may 

be obtained by reference to the following Detailed Description when taken in conjunction with the 
accompanying Drawings wherein: 

FIGURE 1 is an exemplary systemblockdiagramforperformingtheprinciplesofthepresent 
invention; 

10 FIGURE 2 is an exemplary flow diagram describing operations of a multi-dimensional 

kinematic and dynamic modeling simulation having tolerance variations according to the principles of 

the present invention executed by the system of FIGURE 1; 

FIGURES 3 a and 3b (collectively FIGURE 3) is a detailed flow diagram for performing the 

operations according to FIGURE 2; and 
1 5 FIGURE 4 is an exemplary graphical user interface (GUI) for building, displaying, and 

simulating a multi-dimensional kinematic model according to FIGURE 3. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

The principles of the present invention will now be described more fiilly hereinafter with 
reference to the accompanying drawings, in which exemplary embodiments of the present invention 
are shown. This invention may, however, be embodied in many diflFerent forms and should not be 
5 construed as limited to the embodiments set forth herein; rather, these embodiments are provided so 
that this disclosure will be thorough and complete, and will My convey the scope of the invention to 
those skilled in the art. 

Multi-dimensional kinematic and dynamic modeling is an important aspect to mechanical 
engineering design and development projects. Kinematic and dynamic modeling provides a design 

1 0 engineer invaluable information for designing a robust mechanism. By modeling a mechanism and 
performing a dynamic simulation, the design engineer may determine structural weaknesses, accuracy 
of specification, and tolerance robustness, for example. A dynamic modeling system that provides 
results for the design engineer to comprehensively determine performance of the mechanism may 
include a process for simultaneously simulating kinematic and tolerance variation analysis. The 

1 5 simulation may include using a Monte Carlo random variable technique in varying either or both of 
the kinematics and tolerances. 
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FIGURE 1 is an exemplary system block diagram 1 00 for performing the principles of the 
present invention. A computing system 1 05a includes a processor 1 1 0 coupled to a memory 115. 
The processor 1 10 is further coupled to an input/output (I/O) device 120, which may further be 
coupled to a storage medium 125. The processor 1 10 may be a single processor or multiple 
5 processors operating in parallel A control device 127, such as a computer mouse, keyboard, and/or 
% touchpad, may be coupled to the computing system 105 a for controlling a software program 129 

operating on the processor 1 1 0. A display 1 3 0 is connected to the I/O device 1 20 in the computing 
^ system 105a for displaying information generated by the processor. 

r Inoperation, according to the principles ofthe present invention, the computing system 105a 

H 1 0 is utilized to generate and display a kinematic model of a mechanical component 1 3 5 . To generate 
the kinematic model of the mechanical component 1 3 5 , the processor 1 1 0 executes the software 
program 129 stored in the memory 1 15 as loaded from the storage medium 125. A user of the 
computing system 105a controls the execution of the software program 129 by using the control 
device 127. The kinematic model ofthe mechanical component 135 may be generated by selecting 
15 and grouping, via the control device 127, piece parts as rendered on the display 130. Themodelmay 
define orientation and interaction of the mechanical component 135. Alternatively, the mechanical 
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component 135 may be generated by othermeans as knownin the art, suchasusing a keyboard and 
alphanumeric expressions (e.g., a software program). 

The computing system lOSamay be coupled to anetwork 140 whereby data packets 145 
may be communicated across the network 140 to aremote computing system 105b. Thenetwork 
5 140 may be a local area network (LAN), wide area network (WAN), or the Internet, for example. 
The computing system 105a may be a server for executing the software program 1 29 as accessed 
and controlled by the remote computing system 105b. 

FIGURE 2 is an exemplary flow diagram describing operations of a multi-dimensional 
kinematic and dynamic modeling simulation according to the principles of the present invention 
1 0 executed by the system 1 00 of FIGURE 1 . The process starts at step 205 . At step 2 1 0, a kinematic 
model is defined to form a mechanical mechanism (e.g., automobile engine valve) . In defining the 
kinematic model, piece part components may be arranged or connected using predefined software 
tools, such as graphical or programming tools. For example, a graphical user interface (GUI) may 
have predefined piece part tools for rods, gears, hinges, etc., and the piece parts maybe coupled in 
15 a graphical manner. 

At step 215, reference frames (e.g., multi-dimensional coordinate systems) may be applied 
to the kinematic model . The reference frames generally include an origin or base frame that operates 
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as an absolute reference and reference frames at joints to operated as relative references. Similar to 
the defining of the kinematic model at step 2 1 0, graphical or programming tools may be used to apply 
the reference frames to the kinematic model. The reference frames are generally applied at 
connection points (i.e., joints) of the piece parts that are used to form the kinematic model. By 
5 defining the reference frames, a simulator may be utilized to operate the kinematic model having a 
knowledge of the relative and absolute positions of the piece parts. 

At step 220, tolerance parameters are applied to the kinematic model. The tolerance 
parameters generally take the form of +/- X, where X is a decimal number. A technique known in 
the art as geometric dimensioning and tolerancing (GD&T) is a generic type of generating tolerance 

1 0 parameters and may be used to form the tolerances. The tolerance parameters may be aligned with 
the reference frames located at the joints of the piece parts. Alternatively, the tolerance parameters 
may be defined with different reference frames having different origins than the reference frames. The 
order of applying the tolerances (e.g., tolerance parameters before frames or frames before 
tolerances) does not affect results of the simulation. It should be understood that the tolerance 

1 5 parameters generally have the same number of multi-dimensions as the kinematic model, so that if the 
kinematic model is three-dimensional, the tolerance parameters are defined in three dimensions, 
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At St ep 22 5 , the kinematic model including the tolerance parameters may be dynamically and 
simultaneously simulated. In performing the simulation, a Monte Carlo random model simulation 
technique may be employed. The Monte Carlo simulation technique generally applies a random or 
pseudo-random function to a selected variable. For example, the simulation may select one tolerance 
5 parameter, in one or more dimension, to apply a pseudo-random function having a Gaussian or other 
distribution. 

During the simulation, the simulation drives the kinematic model in a predetermined movement 
(e.g. , engine valve moving up and down) while also applying the pseudo-random function to the 
tolerance parameter. The simulation may drive the kinematic model through the same movement 

1 0 while applying the pseudo-random function to different tolerance parameters on an individual basis. 
Alternatively, multiple tolerance parameters may be appUed a pseudo-random distribution during the 
simulation. It should be understood that other variations of the simulation may be employed to 
produce different statistical effects. However, by performing a dynamic simulation that simultaneously 
varies both the kinematic model and the tolerance parameters, a design engineer may view and obtain 

1 5 important information of the mechanical mechanism that is not po ssible by dynamically varying the 
kinematic model and tolerance parameters separately (i.e., serially). The process ends at step 23 0. 
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FIGURES 3aand 3b (collectively FIGURE 3) is a detailed flow diagram 300 forperforming 
the operations according to FIGURE 2. The process starts at step 3 02. A number of general steps 
are shown in dashed boxes. At step 304, individual piece parts may be received as designed by a 
design engineer using a computer aided design (CAD) software program. Alternatively, the principles 
5 of the present invention may be incorporated into a CAD system. At step 3 06, reference frames, 
including joint reference frames, coordinate systems, and datum reference frames, maybe applied 
to the mechanical mechanism may be applied by selecting datum points on the piece parts. 
Additionally, tolerance parameters may be applied to the mechanical mechanism by applying the 
tolerance parameters to the piece parts. Using the reference frames, the piece parts may be 
1 0 assembled into a final mechanism at step 3 08 . At step 3 1 0, a mechanical mechanism is created from 
the piece parts. 

As understood in the art, the number of joint constraints for each joint depends on the joint 
type. For example, aweldjoint has sixjoint constraints. Revolution, roller, cam or prismatic joints 
have five joint constraints. Cylindrical or universal joints has four joint constraints. Spherical or 
1 5 planar joints have three joint constraints. Finally, a bearing joint has two joint constraints. 

At step 312, the simulation is initialized and commenced. The initialization may include 
initializing a Jacobian matrix (Jx) based on the number of j oint constraints, drivers, and piece parts 
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involved. Drivers represent a power source that determines the motion of the mechanism and are 
added to joint axes to set the mechanism in motion in a specific way. At step 3 1 4, the Monte Carlo 
simulationis executed. The simulation may include performing a Monte Carlo simulation for individual 
or multiple tolerance variables simultaneously while performing a dynamic movement of the kinematic 
5 model for the mechanical mechanism. The simulation may end upon fault condition being met, the 
final position has been reached, or the simulation duration pre-specified has ended. 

Specifically, the methodology of the present invention is shown with regard to the particular 
steps shown in solid blocks. In general step 3 04, the component information is read into separate 
part windows at step 3 1 6 fi-om the CAD software. The component information may be imported in 

1 0 the form of International Graphics Exchange Syntax (IGES) files into a operating GUI. At step 3 1 8, 
materials for all piece parts may be defined to determine the properties of parts fi-om a database that 
store the piece part data. At step 320, mass properties including center of mass may be computed. 
For each component in the assembly, the mass properties computed are mass, volume, surface area, 
and moment of inertia, and center of mass, for example. At step 322, geometrical features for 

1 5 emulation are selected fi*om CAD or part model points. Specifically, the geometrical features on the 
piece part may be emulated by either picking vertices on CAD data, by entering the coordinates of 
the model points, or by selecting features in the CAD data. 
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In general step 3 06, reference frames or coordinate systems are then constructed from the 
selected datum points at step 324, and valid datum reference frames designated as joint reference 
frames by selecting the model points created in step 322, if applicable, at step 326. Where 
appropriate, datum reference frames created in step324 may also be used as joint reference frames. 
5 At step 328, deviations and tolerances on respective parts are then defined along with 

additional geometry, e.g. , vectors, needed to define tolerances and measurements at step 3 3 0, To 
simulate variation of a feature, additional geometry required, such as a line, vector, or axis is also 
created at step 330. Further, to define a measurement, additional geometry required, such as a line, 
vector, or axis may also be created at step 330. In general step 308, assembly coordinate systems 

10 are then assigned on the parts or subassemblies at step 332. The datum reference frames or 
coordinate systems defined in step 322 may be designated as assembly coordinate systems in the 
parts or subassemblies. The datum reference frames or coordinate systems may also be used to 
designate the joint type between the parts/sub-assemblies resulting in an assembly. 

At step 334, a determination may be made whether any assembly tolerances are to be 

15 applied. If yes (i.e., assembly tolerances are to be applied), then deviations or tolerance on 
components are defined with respect to datum reference frames or coordinate systems at step 336; 
otherwise, where applicable, drivers are defined at step 3 3 8 . At step 340, magnitude, point or area 
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of application, and direction of force applied at the assembly may be defined. At step 342, 
measurements on subassemblies and final assembly are defined to evaluate pre-determined build 
objectives. 

At step 344, a determination is made whether an initial assembly configuration exists, i.e. , 
5 whether the initial assembly configuration is possible based upon the information given. If not, then 
a further determination is made as to whether the assembly definition requires refinement at step 346. 
If not, then the simulation ends at step 348; otherwise control reverts back to step 324 for fiirther 
datum and reference fi*ame selection. 

In general step 310, if the assembly configuration exists, then a virtual assembly is generated, 
10 e.g., using a Monte Carlo simulation at step 350. The virtual assembly may not satisfy all joint 
constraints for a closed-loop system due to variations. The general coordinates of the parts in the 
assembly may then be used for an initial trial (i.e., x^^^ ) for kinematics analysis. Details regarding the 
general coordinates are discussed afl;er discussion of FIGURE 3B. At step 352, a Jacobian matrix 
defining relative positions of the piece parts is initialized. The size of the Jacobian matrix is initialized 
1 5 based on the number of joint constraints, drivers, and piece parts involved. Again, the number of joint 
constraints in each joint depends on its type. 
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In general step 3 1 2, the kinematic and dynamic simulation analysis starts at step 3 54 having 
a given time (driver positions) for the assembly configuration generated by the Monte Carlo 
simulation. For a new assembly configuration, time (driver positions) are reinitialized. At step356, 
static analyses are performed to determine the assembly float in the joints. Results depend on the 
direction of reaction in the joints. 

In general step 3 1 4, requisite joint constraints are added to the Jacobian matrix at step 358: 
For each j oint, joint constraints may be associated with proper elements of the Jacobian matrix J(x) 
and position vector F(x) based on the joint type and its associated parts. It should be noted that the 
column of the Jacobian matrix where the j oint constraints to be applied depends on its connected 
parts. At step 360, driver constraints are added to the Jacobian matrix: . For each driver, driver 
constraints are placed to proper elements ofthe Jacobian matrix J(x) and position vector F(x) based 
on the driver type and its associated part s. It should be noted that the column ofthe Jacobian matrix 
where the driver constraints to be applied depends on its related parts. 

At step 3 62, flexibility constraints of part s are added to the Jacobian matrix. For each part 
(except ground), flexibility constraints are placed to proper elements ofthe Jacobian matrix J(x) and 
position vector F(x). It should be noted that the column ofthe Jacobian matrix where the flexibility 
constraints to be apphed also depends on the part itself At step 3 64, step changes for the general 
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coordinates are computed, for example, usingNewton-Raphson and Gaussian elimination methods. 
The step changes of general coordinates Dx^^^ are computed by solving a set of linear equations (A6), 
presented below, using Gaussian elimination and backward substitution with total pivoting. The 
advantage of total pivoting allows for taking singularity into account, such as multiple solutions. 
5 At step 3 66, after updating the general coordinates using equation (A7), a determination is 

made at step 368 as to whether the desired precision, allowable constraint tolerances, and flexibility 
of piece parts are met. For example, the step changes of general coordinates Dx^^^ are checked to 
determine whether desired precision is met. Also, at step 3 68, the constraint equations are checked 
to determine whether allowable tolerances and conditions ofthe flexibility of piece parts are satisfied, 

1 0 If any condition is not met, control reverts back to step 3 5 8 for re-evaluating the Jacobian 

matrix J(x) and position vector F(x) using the updated general coordinates (i.e., performed Newton- 
Raphson' s iterations). However, if allowable iterations are reached and the conditions are still not 
met, it means that no solutions exist for the current assembly configuration and driver positions. And, 
general coordinates of piece parts are restored to the previous convergent solution. Otherwise, 

1 5 if all conditions described above are met at step 368, a further determination is made at step 370 as 
to whether the kinematic/dynamic simulation is to be finished. 
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If it is determined at step 3 70 that the simulation is not to be finished, then the time step is 
advanced, and the kinematics and dynamics simulation analysis is repeated for the current assembly 
configuration at step 372 and control reverts back to step 354; otherwise, a further determination 
is made at step 3 74 as to whether the Monte Carol simulation is finished. If not, control reverts back 
5 to step 350; otherwise, the simulation ends at step 348. 

According to Nevrton-Raphson's Method for nonlinear systems, a given set of nonlinear 
equations may be defined: 

i (Al) 

The solution may be converged by the iteration procedure evolved fi*om initial solution x^^^ 
10 as follows: 

^(k) ^ ^(k^i) _j(^(k~j) J (A2) 
where A: > 1 and the matrix J(x) is called Jacobian matrix, defined as: 
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J(x) = 




(A3) 



and vectors x and F(x) are defined as follows: 



7,W 
/2W 



F(x) = 



(A4, A5) 



Therefore, a set of nonlinear equations will be reduced to solve a sequential of linear system 
equations. 
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J(x^''-'> )Axf'> = -F(x''''> ) (A6) 

where 

AX<'> = X^'^-X^'-'^ (A7) 

Gaussian elimination and backward substitution with maximal or total pivoting is then applied. 
5 Gaussian elimination and backward substitution may be performed by the following 

procedure; 



Given a set of linear equations 



Ej :a,jx, +0,2X2 +--- + a,„x„^bj 

E2 : a^jXj + a,,x, + • ■ • + a,„x„ =6, 
: : : : : 



An augmented matrix j? n then formed: 



10 



a = [a b] 



^21 ^22 



«2« ^2,^7 



(B2) 
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where A denotes the matrix formed by the coefficients and the entries in the («+ 1 )-th column are 
the values of b that is a,,„+i = b, for each i=l,2, ...,n. Provided that Ou ^ 0, the operations 
corresponding to . - (a^, /a,, )Ej ) ^ J are performed for each j = 2,3,...,n to eliminate 
the coefficient ofxi in each of these rows. Follow a sequential procedure for/ = 2, 3, . and 
5 perform the operation (e^ -(aJa„)E,)-^(E^)foTeachj = i+l,i+2, provided that a„ = 

0. Thiseliminates the coefficient ofx,in£, for all valuesof/= 1,2, (2<j<n). Theresulting 
augmented matrix has the form: 



A = 



0 «22 - «2„ 



0 «„.„^7 



(B3) 



where, except in the first row, the values of a,^ are not expected to agree with those in the original 

10 matrix A . This matrix represents a linear system with the same solution set as system (Bl). Since 
the equivalent linear system is triangular: 



(B4) 
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backward substitution can be performed. Solving the n-th equation for x„ gives 



(B5) 



Solving the (w-l)-th equation for x„.i and using x„ yields 



^n-, = (B6) 



5 and continuing this process, all variables can be obtained by 



X. = 



a.. 



(B7) 



for e^chi =11-1,11-2,'" J J 



However, the procedure will not work if one of diagonal elements a^^ is zero during elimination 
procedure. Therefore, pivoting strategies are required. 
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Maximal (or total) pivoting maybe performed by using the following procedure. Maximal 
pivoting at the k-th step searches all the entries ajj, for i = k, k+1 , . . ., n, and j = k, k+1 , . . ., n to find 
the entry with the largest magnitude. B oth row and column interchanges are performed to bring this 
entry to the pivot position. 

In further discussion of the Euler parameters approach (as utilized in step 366), general 
coordinates are used for relative position of any piece part with respect to a global coordinate system 
of the whole assembly (or ground - non-moving part) is represented by transformation matrix as 
follows: 

"4f+e|-0.5) 2{e,e,-e,e,) 2{e,e, + e,e,) tj 
2{e,e,+e^e,) 2{ef + e', -O.5) 2{e,e,-e,e,) 
2{e,e,-e,e,) 2{e,e,+e,e,) 2{e^ + e', - O.5) 

0 0 0 1 



10 where rotation matrix are composed of e^, 62, e^, called Euler parameters and ty, t, are 
translation components. These seven variables are denoted as general coordinates of the piece part. 
Here, an additional variable is introduced in rotation matrix. If the part is a rigid body, additional 

constraint, e/ + + + = i , is applied. 
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For agivenjoint Joint constraint equations, which are expressed by general coordinates of 
two connected parts, may be derived by relating two joint reference frames that are attached to the 
connected parts, respectively. For example, a revolution] oint contains five equations to constrain 
three translation degrees of freedom (OjOj - 0) and two rotation degrees of freedom (e,j 1 e,i & 
Cjj 1 Aprismaticjoint contains five equations to constrain three rotation degrees of freedom (e^j 
1 e^i & Cyj 1 e,i & e,,; 1 Cyj) and two translation degrees of freedom (e^i 1 OjOj & Cyj 1 OiOj). 
Although types of joints may be different, constraint equations are similar. Therefore, a set of 
nonlinear joint constraint equations expressed by general coordinates of related parts may be easily 
obtained and so are their Jacobian components. 

In general, there are two types of driver constraint equations. One is for rotation driver (Z 
(Cxi, e^j) and the other is for translation driver (OjOj • ej. Similarly, they can be easily derived by 
relating two joint reference frames that are attached to each connected part respectively. By using 
the Euler parameters approach, the following may be achieved: 

1 . Number of independent variables (general coordinates) and constraint equations 
depends only on the number of joint constraints, drivers, and parts involved. The 
independent variables are independent of the number ofloops existed in the assembly 
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and how joints located in a loop as in the vector loop approach. Therefore, the Euler 
parameters approach may be appUed for any generic mechanisms. 

2 . Flexibility (i.e. , the reverse of stiflBiess) of a piece part may be taken into account by 
introducing Euler parameters. The more flexible the piece part, the larger an offset 

5 value is between the constraint, + e-i + e-^ + and 1. It is impossible to 

consider flexibility of a piece part in the vector loop approach. 

3 . The ways that joint constraint equations are set up allow a user to control joint 
mismatch or axis misalignment. This feature is used when integrating kinematics 
analysis with dimension tolerance variation simulation. It is impossible to implement 

10 this feature in the vector loop approach. 

It should be understood that upon adopting the Euler parameters approach, a user may 
perform the following modifications to the process ofFIGURE 3 to take into account special cases 
to improve robustness of the system. 

1 . Use an inverse calculation of the Euler parameters from a given transformation matrix to 
15 take into account the possibility of a special singular case. 
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2. Take into consideration the relationship between a local transformation matrix and a 
global transformation matrix for a piece part into considerations so that more generic assembled 
sequences may be applied. 

3 . Consider equivalent angular solutions during iterations. 

5 4. Modifythe Gaussian elimination and backward substitution with total pivotingtotakeinto 

account a case that may have multiple solutions. 

5 . Allow a user to have control over the piece part flexibility and the constraint tolerances. 
FIGURE4 is an exemplary graphical user interface (GUI) 400 for defining, displaying, and 
simulating a multi-dimensional kinematic model 405 according to HGURE 3 . The kinematic model 

10 40 5, which may be imported fi-om a CAD system, is shown in window 410. Elements 4 1 5a-4 1 5h 
in window 420 may be utilized to define joint types for placing reference frames on the joints 
connecting the piece parts of the kinematic model 405. Window 425 provides for setting and 
adjusting parameters for controlling the simulation. As shown, the parameters may include: (i) start 
time, (ii) duration, (iii) increment, and (iv) analysis precision. Other parameters may also beincluded. 

1 5 In operation, the design engineer may execute the simulation and view and/or obtain dataresulting 
from dynamically moving the kinematic model 405 while simultaneously performing tolerance 
variations. 
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The previous description is of apreferred embodiment for implementing the invention, and 
the scope of the invention should not necessarily be limited by this description. The scope of the 
present invention is instead defined by the following claims. 
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